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Rigid cluster decomposition reveals criticality in frictional jamming 
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We study the nature of the frictional jamming transition within the framework of rigidity percolation theory. 
Slowly sheared frictional packings are decomposed into rigid clusters and floppy regions with a generalization 
of the pebble game including frictional contacts. We discover a second-order transition controlled by the emer¬ 
gence of a system-spanning rigid cluster accompanied by a critical cluster size distribution. Rigid clusters also 
correlate with common measures of rigidity. We contrast this result with frictionless jamming, where the rigid 
cluster size distribution is noncritical. 


The interplay of constraints, forces, and driving gives rise 
to the jamming transition in granular media. It is now well- 
established that the frictionless jamming transition has char¬ 
acteristics of both second- and first-order transitions. Both 
the average coordination number and the largest rigid clus¬ 
ter size jump at the transition, yet there exists a diverging 
lengthscale [1—4], Frictional jamming is more puzzling: The 
hysteresis observed in the stress-strain rate curves of stress- 
controlled flow simulations [5-8] and experiments [9] has lead 
to an interpretation as a first-order transition. Yet, signs of 
second-order criticality appear when treating the fraction of 
contacts at the Coulomb threshold as an additional parame¬ 
ter [10-12], 

To elucidate the frictional jamming transition from a micro¬ 
scopic viewpoint, we extend concepts and tools from rigidity 
percolation, i.e., the onset of mechanical rigidity in disordered 
spring networks [13-16], to frictional packings. The former 
is driven by the emergence of a system-spanning rigid cluster 
that can be mapped out (in 2d) using the pebble game [17], an 
improved constraint counting method that goes beyond mean- 
field by identifying redundant constraints. We, for the first 
time, implement a generalized pebble game for 2d frictional 
systems and use it to identify rigid clusters in very slowly 
sheared packings. As we show below, this allows us to iden¬ 
tify a second-order rigidity transition and to link stresses and 
nonaffine motions to the microscopic stmcture of frictionally 
jammed packings. 

Generalized isostaticity: To establish context, we first re¬ 
view the application of Maxwell constraint counting to jam¬ 
ming [18]. For A particles in d dimensions and a mean num¬ 
ber of contacts per particle z, interparticle forces yield dzN /2 
constraints. Since each particle has \d(d + 1) translational 
and rotational degrees of freedom, there are y(N — l)d(d + 1) 
total degrees of freedom (subtracting out global degrees of 
freedom). When these match the force constraints, we arrive 
at the isostatic criterion, or dzN /2 = \{N — 1 )d(d + 1). In 
the limit N —> °°, 7 iso = d+ 1 for frictional granular materials. 
For frictionless packings, we ignore rotations and obtain the 
familiar Ziso = 2 d. 

Despite being mean field, i.e. neglecting spatial correla¬ 



FIG. 1. The pebble game and rigidity percolation for a sample 
N = 1024 frictional packing, a - Force chains, with anisotropic 
stresses due to Lees-Edwards boundary conditions, b - Contact net¬ 
work, with frictional (black) and sliding contacts (red), c - Result 
of the (3,3) pebble game where frictional and sliding contacts have 
been mapped to double (thick) and single (thin) bonds, respectively. 
Red are pebble-covered bonds, green are overconstrained bonds, and 
colored disks are leftover pebbles, d - Rigid cluster decomposition: 
This packing is partially rigid and consists of three rigid clusters 
(black, purple and green) and regions of floppy bonds (gray). 

tions, isostaticity works seemingly well to locate the jamming 
transition in static frictionless systems [1]. For frictional sys¬ 
tems, however, numerical and experimental evidence point to 
a range d + 1 < z < 2d at the transition, with a matching den¬ 
sity range from random loose packing [19] to random close 
packing. To resolve this conundrum, a generalized isostaticity 
criterion was introduced [10, 11], that accounts for contacts at 
the Coulomb friction threshold providing one less constraint: 

z? so = (d+l)+2n m /d, (1) 

where n m is the mean number of such fully mobilized con¬ 
tacts per particle. Equation (1) describes a line of transition 
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points interpolating from z = d+ 1 at n m = 0 , corresponding 
to the friction coefficient p = °° limit, to n m = 1 at z; so = 4, 
corresponding to the p = 0 limit [ 10 ]. 

Simulation. To obtain frictional packings near jamming, 
energy-minimization is not an option for the non-conservative 
frictional interaction. Instead, we implement a common ex¬ 
perimental and simulation protocol, simple shear at strain rate 
7 in Lees-Edwards boundary conditions. As we decrease the 
strain rate, we expect to reach a “quasistatic” limit, where the 
system transitions through a series of (nearly) force-balanced 
states so that the mechanical properties of these states are rel¬ 
evant, and we can use the tools of rigidity percolation. 

We simulate systems of N = 64 — 4096 polydisperse parti¬ 
cles in two dimensions interacting according to the Cundall- 
Strack law [20]. In addition to harmonic purely repulsive 
normal forces /„ = k„8, (8 is the particle overlap), it com¬ 
prises an incremental tangential force df lm = k t dt, where dt 
is the amount of tangential sliding since the establishment of 
the contact. The friction force is constrained by the Coulomb 
criterion | f, \ < pf n ; once the threshold is reached, it continues 
to slide at \f = fif, until the direction reverses. The energy 
dumped into the system by shearing is dissipated through lin¬ 
ear viscous damping forces, ffj = — £(v,- — v ; ), and a small 
amount of rotational individual damping. We work in scaled 
units with mean particle radius (r) = 1 and unit stiffness 
k n = k t = 1. Most of the results below are in the low-friction 
p = 0.1, low-damping £ = 0.1 limit, for A = 1024 particles 
except where specified otherwise. We report additional re¬ 
sults for high p = 1 and high £ = 1 in Figures S4 and S5, 
respectively. Systems are strained for T = 10 6 time units in 
all cases, corresponding to a strain of 10 system lengths for 
7 = 10~ 5 strain rate, and 1 system length for 7 = 10~ 6 . 

Rigidity percolation. We decompose packings into rigid 
clusters using a (k = 3,l = 3) pebble game [17]: First, we as¬ 
sociate a pebble with each of the k = 3 degrees of freedom of a 
particle. We then build a constraint network based on the con¬ 
tact network where a contact imposing n constraints translates 
to n bonds, i.e., fully mobilized contacts with \f t \ = pf n yield 
one bond while contacts with \f,\ < pf„ yield two (Fig. SI). 
We then recursively move these pebbles along bonds in the 
constraint network to assign pebbles to bonds, where each 
pebble-covered bond is equivalent to a degree of freedom be¬ 
ing independently constrained. We do so until l = 3 or more 
pebbles cannot be assigned. Finally, we map out the rigid 
clusters for each contact network. This algorithm is an exten¬ 
sion of the (2,3) pebble game used in the frictionless case. 
Please see the SI and Fig. S2 for more details. 

In Fig. 1, we show the four stages of moving from the sim¬ 
ulated packing to the rigid cluster decomposition. The force 
chains in Fig. la correspond to the constraint network in Fig. 
lb and lc, with frictional contacts (double bonds) in black, 
and sliding contacts (single bonds) in red. This network forms 
the basis for the (3,3) pebble game in Fig.lc. In our exam¬ 
ple, it is possible to assign pebbles to most of the bonds in 
the network without overconstraining it. In agreement with 
the beyond-mean field nature of Faman’s theorem (see SI), a 


number of leftover pebbles also remain, especially at rattlers 
buth also at connected particles (colored circles). Finally, in 
Fig. Id, we decompose the system into rigid clusters, that is 
connected regions where no more than the 3 pebbles linked 
to the global translations and rotations can be found. We find 
three rigid clusters including a large, system-spanning one. 
The remaining bonds are floppy, i.e. not rigid with respect 
to any of their neighbouring bonds. Our example system is 
globally rigid, in spite of an average z. below the generalized 
isostaticity criterion. We show that such configurations are 
generic below. 

Results: We first address how the global properties of the 
system, including mean stresses and distance from isostaticity, 
depend on density, strain rate and system size (Fig. 2). Exper¬ 
imental [9] and simulated frictional systems [5-8] report a 
hysteresis loop in the stress-strain relations, through a proto¬ 
col that includes either a strain rate ramp or constant stress 
driving. As a function of density, (j) = 0.8225 — 0.84, we see 
a transition through jamming, as evidenced by the pressure 
distribution shifting from a peak at 10 4 (in units of over¬ 
lap) to a peak at 10 ~ 2 (Fig. 2a; in Fig. S3 we show that p 
and a xy are equivalent). Since we perform a constant strain 
rate simulation, we do not observe hysteresis. We instead 
find intermittent flips between jammed and unjammed states 
for the intermediate density samples, in particular 0 = 0.8275 
and <j> = 0.83 (see Fig. 3). When lowering strain rates, a gap 
in pressure opens between low and high densities, consistent 
with an approximately Bagnold scaling p ~ y 1 2 dominated 
by viscous damping forces below jamming, and the appear¬ 
ance of a yield stress above jamming (Fig. 2b). 

We focus our attention on the tuning parameter dz = z — z!" n . 
i.e. the distance from generalized isostaticity. For each strain 
rate, all of the data for different packing fractions and system 
sizes collapses onto a unique curve p(dz) (Fig. 2c). For the 
two lower strain rates, we see two regimes - a rapid drop be¬ 
low dz~ 0 which depends on 7 , separated from a more grad¬ 
ual, universal increase at dz > 0. The small, positive values 
of the pressure for dz < 0 strongly depend on £ (see Fig. 
S5), indicating again that viscous damping forces dominate 
this regime, consistent with recent results for shear thicken¬ 
ing in suspensions [21]. When we visualize our system in 
a two dimensional n m — z. phase diagram (Fig. 2d), jammed 
(blue) states defined here by p > 10 2 exist predominantly 
in the stable region of the phase diagram below the stability 
line, while unjammed (red) states mostly exist in the unsta¬ 
ble region. There are however significant fluctuations in the 
transition region with some jammed states below the stability 
line, suggesting that this mean-field criterion is insufficient 
(Fig. 2e). System trajectories have roughly equal fluctuations 
in z and n m . unlike the avalanching system of [ 12 ], where fluc¬ 
tuations along n m were more prominent, showing that the na¬ 
ture of the trajectories depends on the type of driving. 

We now present the results of the rigid cluster decomposi¬ 
tion using the (3,3) pebble game. To demonstrate the struc¬ 
tural importance of the rigid clusters, we first correlate the 
time series of rigidity and pressure in the region where we ob- 
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FIG. 2. a - Pressure distributions for 8 densities (j) = 0.825 — 0.84 accross the frictional jamming transition, for N = 1024 and y = 10~ 5 . 
b - Stress-strain relations for the same densities and N, and strain rates y = 10~ 5 — 10~ 7 . c - Correlation between distance from generalized 
isostaticity dz = z— (3 + n m ) and pressure (log-scale), d - Generalized isostaticity: Probability density of the system in ( n m ,Zi S0 ) phase space, 
combined from densities <p = 0.825 — 0.845; N = 1024 and y = 10~ 6 . The red area has p < 10 -3 and is unjammed, while the blue region has 
p > 10~ 3 and is jammed, e - Fluctuations for the same data set, each dot is a jammed (blue) or an unjammed (red) packing. 



FIG. 3. Temporal correlation between rigidity and stresses. Left: 
Pressure trajectories for four densities close to jamming (log-scale), 
followed by the fraction of the system that is rigid, (frac), the x-length 
of the largest cluster and the y-length of the largest cluster. Right: 
Time correlation between pressure and rigidity (frac(f) log p(t+dt)). 

serve intermittent behaviour (see Fig. 3). In the second, third 
and fourth panel, respectively, we show the fraction frac of the 
bonds belonging to a rigid cluster, the x-extent L x and the y- 
extent L y of the largest cluster (normalized by system size L). 
Rigid systems are characterized by frac sa 1, and system span¬ 
ning clusters in the x and y directions. All three measures cor¬ 
relate with pressure and with each other, though switches be¬ 
tween globally rigid and floppy states are significantly faster 
than pressure changes. On the right (Fig. 3b), the correlation 
function (frac(f)logp(f +dt)} shows strong, symmetric corre¬ 
lations for the mostly unjammed runs (j) = 0.8225 — 0.83, and 
a slight asymmetry indicating that pressure follows rigidity, 
for the jammed runs 0 = 0.8325 — 0.84. 

Fig. 2 demonstrates that dz is an appropriate mean-field 
parameter. In light of Fig. 3, we conclude that the spatial de¬ 
composition indeed plays a role in the macroscopic response 
of the system. We now ask how rigid cluster analysis can help 
uncover the nature of the frictional jamming transition. In Fig. 
4, top row left, we show the rigid cluster size distribution p(n), 
where n is the number of bonds in a rigid cluster, across the 
frictional jamming transition. We observe a change character¬ 
istic of a second-order phase transition, with a cluster size dis¬ 


tribution which broadens approaching the transition, and then 
the emergence of a system-size percolating cluster above the 
transition. At the transition, we observe a power-law distribu¬ 
tion with an exponent a ~ —2.5. For comparison, for connec¬ 
tivity percolation in two dimensions a = —187/81 = —2.31 
[22] and a self-organized rigidity percolation model yields 
a = —1.94 [23]. To help pinpoint the location of the tran¬ 
sition, we plot the maximum cluster length (L 2 + L 2 ) 1 / 2 /y/2L 
against dz. It approaches unity near dz = —0.15 rather than 
dz = 0, consistent with the picture of rigid and floppy regions 
coexisting in an overall rigid system emerging from Fig. 1. 
So does an equivalent measure, the spanning probability (Fig. 
S7). Moreover, this downward shift survives in both the large 
N (Fig. 4 top right) and the y—> 0 (Fig. S6) limit. 

As comparison, we simulate a frictionless system across its 
frictionless jamming transition [24] with the same protocol as 
our frictional runs (Fig. 4, bottom row). Rotations are irrel¬ 
evant to frictionless disks, so we use a standard (2,3) pebble 
game here (see SI). In stark contrast to the second order tran¬ 
sition discussed above, here we observe the hallmarks of a 
first-order transition: The rigid cluster size distribution is ei¬ 
ther rapidly decaying at low pressures (most packings have 
no rigid clusters at all), or markedly bimodal without any 
intermediate-sized clusters for an order of magnitude. Also, 
the largest rigid cluster length is gapped (Fig. 4 right, inset), 
and depends strongly on system size. This finding is consis¬ 
tent with recent results for frictionless systems where pack¬ 
ings were found to be either fully rigid or fully floppy [2]. 
Again, we observe a downward shift in the transition point, 
though it is not as clear if it survives the large N limit or the 
y —> 0 limit. Recent work demonstrates that the frictionless 
transition point can in fact change with protocol [25]. 

To complete our analysis, we correlate rigid clusters with 
two other microscopic measures of rigidity for our frictional 
packings (see Fig. 5). First we measure the mean bond nor¬ 
mal force scaled by the packing pressure /„/p over bonds be¬ 
longing to either a rigid cluster or a floppy region, ensemble- 
averaged over packings with similar dz (left). We see that /„ 
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FIG. 4. Rigid cluster size statistics. Left we show the cluster size 
distribution for pressures accross jamming, while on the right is the 
scaling of the largest cluster size with dz, for 3 system sizes. Top 
row: frictional systems, N = 64 — 1024 , y = 10~ 6 , 5 runs each at 8 
packing fractions with 160000 packings in the transition region. Bot¬ 
tom row: frictionless systems, idem with 64000 samples at packing 
fractions <j) = 0.835 — 0.85. 

inside the clusters always closely tracks the pressure, confirm¬ 
ing that clusters indeed always bear stress. In contrast, floppy 
regions markedly depend on dz : Below the transition, they 
bear the same load as rigid regions, but then rapidly lose load 
once the percolating rigid structure forms. Finally, the unsta¬ 
ble regions become isolated rattlers which bear no load. We 
interpret this again as viscous forces dominating the stresses 
below the transition, where no spanning cluster can bear the 
load. The gradual decrease of non-load bearing rattlers above 
jamming is well-known in static systems [26]; we put it into 
a dynamical context here for the first time. Spatial correla¬ 
tions in the forces are also found in shear jamming with the 
transition occuring below the isostatic point [9], The second 
measure are the nonaffine motions, which are known to dra¬ 
matically increase approaching the frictionless [27] and fric¬ 
tional [11] jamming transitions. We measure the relative tan¬ 
gential motion of the centers (| |) and tangential sliding at the 

contacts (|c4mg|)- Let r (/ - = r ; (/) — r j(t) be the vector linking 
two neighbouring particles’ centers, t, ; the tangential unit vec¬ 
tor at the contact, Rj the particle radii, and a, their angles (see 
Fig. SI). Then [11,28] 

d t = Tjj * t/ /■ ^tang = L/'t ij (RjCtj T R jCCj j . (2) 

Normalized by the strain rate y, nonaffine motion is sig¬ 
nalled by values above 0.5. For both rigid and floppy regions, 
above and below the transition, motion is strongly nonaffine 
(Fig. 5, right, note log scale). However, displacements in 
floppy regions are much more nonaffine compared to rigid re¬ 
gions on both sides of the transition, culminating in values 
d t /f> 20 for isolated rattlers in rigid packings. The normal 
displacements remain at d, ,/y« 0.5 throughout (see Fig. S3). 
Links between nonaffine buckling and local rigidity have been 
pointed out previously [29]. We also observe a peak in the to¬ 




FIG. 5. Correlation between rigid clusters and mechanical stability. 
Left: Scaled normal force magnitude /„ /p inside (red) and outside 
(green) the clusters, together with the rigid fraction of the system 
(black). Right: Nonaffine motion (tangential d, and contact slipping 
<ftang) inside and outside the cluster, and total. 

tal nonaffine motion accross the transition. 

Conclusions: In sum, we adapt ideas from rigidity per¬ 
colation to characterize the frictional jamming transition of 
slowly sheared packings beyond mean-field level. We show 
that while generalized isostaticity is a good mean-field crite¬ 
rion, spatial correlations do matter: a packing can be jammed 
below global isostaticity if it contains both floppy regions and 
a system-spanning rigid cluster, resembling spring networks 
in this regard. The emergence of such a cluster appears to be 
first-order in frictionless packings, with a sudden jump from 
microscopic to system-wide clusters, but second-order in the 
frictional case, with a power law distribution of cluster sizes 
at the transition. In particular, partial rigidity is unique to fric¬ 
tional packings. 

The key instrument in obtaining those new results is rigid 
cluster decomposition. It allows us to draw connections be¬ 
tween disordered spring networks, where it is the norm, and 
granular packings [2, 30]. By nature, cluster decomposition 
ignores the detail of contact forces, which may explain why 
it only partially correlates with local stresses: presumably, 
arching is only possible around small floppy regions. It does, 
however, account for the dynamical nature of fully-mobilized 
contacts, thus highlighting their central role in frictional pack¬ 
ings. Applied strain leads to internal rearrangements, moving 
through a phase space of rigid, non-rigid and crucially par¬ 
tially rigid packings due to the second-order nature of the tran¬ 
sition. Incorporating this internal variance of accessible states 
through an internal field, Grob et al. were able to explain 
the phenomenology of hysteresis in frictional jamming [8]. 
Our work begins to provide a microscopic basis for such phe¬ 
nomenology. 
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Supplemental material 


GENERALIZED PEBBLE GAME 

The foundation of extending ideas of constraint-counting 
beyond mean-field in two dimensions is Laman’s theorem [1], 
Specifically, it states that a network with N vertices is gener- 
ically, minimally rigid in two dimensions if and only if it has 
2N — 3 bonds and no subgraph of n vertices has more than 
2n — 3 bonds. Here minimal rigidity signifies that should any 
one contact/bond be deleted from the packing then the system 
breaks up into at least two rigid clusters, i.e. there are no re¬ 
dundant contacts and no floppy modes (prior to the deletion). 
Rigid clusters are defined along contacts/bonds and consist of 
connected bonds that are rigid with respect to each other in 
the sense that there is a nonzero energy cost to moving one 
bond with respect to another. In the language of jamming, the 
existence of a spanning cluster that is minimally rigid trans¬ 
lates to the packing being precisely at the jamming transition 
point. Laman’s theorem is purely combinatorial. The power 
of the theorem is that it goes beyond simple constraint count¬ 
ing, which is just a necessary condition, and additionally ac¬ 
counts for whether or not the constraint is redundant since 
such contacts do not contribute to the emergence of rigidity. 

The (k = 2,1 = 3) pebble game [2] is an efficient algo¬ 
rithm to implement Lamans theorem. It also provides us with 
a method for identifying rigid clusters in the packing. In the 
(2,3) pebble game, k = 2 pebbles, corresponding to the lo¬ 
cal degrees of freedom for each particle, are placed on each 
node/particle of the packing/network [2]. These pebbles can 
be moved along the contacts in the contact network in such a 
way to explore all possible subgraphs. If a contact is indeed 
an independent constraint, then a pebble will be assigned to 
that contact. The pebble game stops when there are l = 3 free 
pebbles (not assigned to any contact), corresponding to the 
I = 3 global degrees of freedom, which do not tell us about 
the nontrivial deformations of the packing. 

While the original work on the pebble game focused on the 
k = 2,1 = 3 case [2], Lee and Streinu addressed the general 
case [3]. Here is a description of the general case based on 
their work [3]. To determine if the contact network is un¬ 
der or over or minimally constrained, the rules are as follows: 
(1) No more than k pebbles are allowed on any nodes in the 
contact network, and (2) a bond is labeled as an independent 
constraint when there are at least l + 1 pebbles in total present 
at the two nodes connected by the contact. Allowable pebble 
moves are: (1) Pebbles may be moved from node to node by 
searching a directed graph such that if a pebble is found then 
the bonds along the directed path are reversed and (2) once a 
pebble has been assigned to a bond, the direction of that bond 
is determined by the node the pebble came from. The game is 
then played until all the bonds have been processed. If there 
are 1 free pebbles left over and no unassigned bonds, then the 
contact network is minimally rigid, or at the jamming transi- 




FIG. SI. (a) Top: Schematic illustrating the contact geometry for 
frictional contacts, the force on particle i due to particle i assuming 
particle j is stationary, and the displacement and rotation of particle 
i. Bottom: Mapping to constraint network, (b) Same as (a) except 
for sliding contacts. 

tion. If there are more than / pebbles left over, the network is 
floppy. If there are unassigned bonds with / pebbles left over, 
then network is overconstrained. 

The pebble game can be extended to compute rigid clusters 
by temporarily pinning / free pebbles and then marking two 
particles in contact as rigid. Neighboring disks in contact are 
then explored to free a pebble. If a free pebble is found, the 
neighboring disks and their associated contacts are flexible. 
If the converse, the associated contacts are rigid. Rigid clus¬ 
ters are defined along contacts/bonds and consist of connected 
contacts that are rigid with respect to each other. This proce¬ 
dure can then be used to map out the rigid clusters. A system 
that is not jammed will not have a spanning rigid cluster, so 
the existence of a spanning rigid cluster determines the state 
of the system (fluid or solid). The existence of a spanning 
minimally rigid cluster determines the transition point. 

In order to identify rigid clusters in frictional packings, we 
first must extend the commonly used (k = 2, l = 3) pebble 
game for frictionless disks, where only normal (central) forces 
exist, to frictional disks, where both normal forces and tangen¬ 
tial forces exist. Since there are now two translational degrees 
of freedom and one rotational degree of freedom for each par¬ 
ticle, k = 3. The constraints can be taken into account as ex¬ 
plained above, i.e. once a pebble is assigned to a contact in the 
pebble game, then that contact is considered an independent 
constraint. However, to account for the additional constraints 
due to tangential forces in the frictional case, we introduce a 
second bond for each frictional contact into what we now call 
the “constraint network”. The pebble game then explores the 
constraint network to see if that additional rotational degree of 
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FIG. S2. The (3,3) pebble game for two sample configurations. Top 
row: Three particles connected by frictional bonds become a rigid 
cluster, since there are a total of 3 pebbles left. Bottom row: Four 
particles in a square configuration are flexible, since there are a total 
of 4 pebbles left. 

freedom can be independently constrained. This second bond 
in the constraint network is only added to frictional contacts 
below the Coulomb criterion, where the normal and tangential 
forces are independent of each other. For fully mobilized con¬ 
tacts (contacts at the Coulomb criterion), the tangential and 
normal forces are no longer independent such that only one 
bond in the constraint network is needed. Since the num¬ 
ber of global degrees of freedom in two dimensions is still 
3, we arrive at a (k =3,1 = 3) pebble game where contacts 
below the Coulomb criterion are denoted as double bonds in 
the constraint network and contacts at the Coulomb criterion 
are denoted as single bonds in the constraint network. While 
the (3,3) pebble game has been used to analyze the rigidity in 
body-bar networks [4], it has not yet been used for frictional 
disks. See Fig. SI. 

As a test of this approach to identifying rigid clusters in 
frictional disk packings, one can consider a few small pack¬ 
ings with three or four disks (see Fig. S2). In the friction¬ 
less case, the smallest nontrivial minimally rigid packing (and 
corresponding constraint network) is a triangle. In the fric¬ 
tional case, if we consider a triangle with each contact below 
the Coloumb criterion, then there are two bonds between each 
node in the constraint network (panel a). After performing 
the (2,3) pebble game, there are 3 free pebbles left and so the 
constraint network is minimally rigid (panel b). Note that this 
result is in keeping with three gears in a triangle not being able 
to turn. And if at least one of the “double bonds” in the con¬ 
straint network is removed from the triangle, then the packing 
is flexible since sliding between two of the disks is possible. 
If instead we consider four frictional disks packed as a square 
with each contact below the Coloumb criterion (panel c), then 
there are four free pebbles left after the (3,3) pebble game is 
played and the packing is flexible just as four gears arranged 
in a square can turn (panel d). 

Finally, we address boundary conditions. For two- 
dimensional systems with periodic boundary conditions. 


i? io 3 in’ io 1 io 3 „ io 5 

FIG. S3. Cluster size distributions for the (3,2) and (2,2) pebble 
games for frictional and frictionless systems, respectively. Left: fric¬ 
tional system at y = 10~ 6 and £ = 1.0. Right: frictionless system, 
same parameters as Figure 4. Despite the limited statistics, the fric¬ 
tional transition still appears continuous, while the frictionless tran¬ 
sition retains a gap in cluster sizes. 

global rotations are forbidden by symmetry and there are only 
two global degrees of freedom (the two translations). Clus¬ 
ters smaller than system size do not suffer this restriction, but 
a spanning cluster can only have two, rather than three, left¬ 
over pebbles. In practice and for our purpose, there is little 
difference between the l = 2 and the 1 = 3 pebble game. Al¬ 
though some results are affected quantitatively, the location of 
the transition point remains the same and our main conclusion 
is unaffected: Frictional systems have a continuous cluster 
size distribution (Figure S3, left), while frictionless systems 
show a distinct gap (Figure S3, right). In other words, the 
precise form of the pebble game does not affect the observed 
order of the transition. 


ADDITIONAL RESULTS 

In this section, we present a series of supplementary re¬ 
sults which complement our message of frictional jamming as 
a rigidity transition, by varying the strain rate y, the friction 
coefficient q and the viscous damping coefficient £, and by 
presenting alternative measures of the transition. 

Stress measures. We begin, however, by confirming that 




FIG. S4. Right: Ratio of other stress components to pressure: normal 
forces, xy, xx and yy components of the stress tensor. There are no 
systematic variations except for /„, where the mean contact number 
has an influence. The bars indicate the standard deviation of the fluc¬ 
tuations, they are not error bars of the mean. Right: Normal part of 
the decomposition of the motion inside and outside the cluster. The 
normal displacements are close to the affine prediction d n /y= 0.5. 
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our use of the pressure as a marker of jamming is consistent 
with other possible definitions, in particular the shear stress 
O xy . Figure S4 shows the ratio of four alternative measures 
of the stresses and the pressure (normal force /„, shear stress 
Oxy, diagonal stress components o xx and o yy ). We have indi¬ 
cated the standard deviation of the fluctuations by error bars. 
Except for /„, where the mean contact number matters, this 
ratio is virtually independent of dz , with a finite and constant 
fluctuation amplitude, indicating that all of these quantities are 
equivalent measures of the distance to the transition. 

Nonaffine motion - normal displacements. Complementing 
Figure 5 in the main text. Figure S4 (right) shows the normal 
components of the local displacments, d'f = dr, 7 - • n ;; -, scaled 
by the strain rate. The values both inside and outside the clus¬ 
ters are close to the affine prediction, d n /y = 0.5. 
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FIG. S5. Properties of a frictional system with friction coefficient 
p = 1. Top left: pressure distributions for 8 values of <j> accross the 
jamming transition. Top right: pressure - dz scaling, showing the 
same collapse as for p =0.1. Bottom left: Two dimensional his¬ 
togram of states, colored according to jammed (blue) and unjammed 
(red); criterion for jamming p > 10~ 3 . Bottom right: Scatter plot of 
the maximum cluster size as a function of dz. 


(panel d). Though our statistics is limited, we can clearly see a 
growing cluster when approaching the transition from below, 
and then a spanning cluster appears somewhat below dz = 0. 



FIG. S6. Results for higher viscous friction values £ = 1.0. Left: 
pressure - dz plots; the differences compared to the low friction case 
^ = 0.1 are at unjammed dz, where there are larger contributions to 
the viscous stresses. Right: Rigidity transition scaling for three dif¬ 
ferent system sizes. The result is essentially unchanged from Figure 
4. 

Influence of £. Figure S6. The viscous damping coefficient 
£ chiefly affects the viscous part of the forces, which is in¬ 
cluded in our stress measurements. Intuitively, the mechanism 
of the transition should not depend on £ for low enough strain 
rate since the viscous forces /j? ~ £(v 7 — v,-) ~ C7 ,; aff g° t° 
zero (the last term is the nonaffine parameter). We find that 
while the value of £ strongly increases the pressure values be¬ 
low jamming, consistent with the viscous forces dominating 
this regime, the pressure in the jammed region is virtually un¬ 
affected (left). The transition curve (right) strongly resembles 
the £ = 0.1 case in Figure 4, albeit somewhat smoothened out. 




Influence of p. In Figure S5, we present equivalent data 
to the main article for p = 1, i.e. in the high friction limit 
as opposed to the p = 0.1 low friction limit. From the gen¬ 
eralized isostaticity point of view, the main influence of the 
higher p is a smaller number of fully mobilized contacts, con¬ 
sistent with [5, 6]. In the n m -z. diagram (panel c), jammed 
configurations have close to no sliding contacts, and rest just 
within the stable region. Unjammed configurations have more 
sliding contacts, and also significantly fewer contacts over¬ 
all. As shown in panel a, the jamming transition shifts to 
lower tj) - from our limited resolution it lies in the range 
0.8 < (f>j < 0.805. Unlike for p = 0.1 in Figure 3, we do 
not observe runs which flip between jammed and unjammed 
states, possibly since we are not close enough jamming. If we 
plot pressure vs. dz, our data collapse on nearly an identical 
curve to the equivalent N = 1024, y = 10 6 line in Figure 2. 
Finally, we have repeated the rigid cluster analysis for p = 1 


FIG. S7. Rigidity transition properties of the frictional system as 
a function of strain rate: for decreasing strain rates, the transition 
becomes sharper. Left: length of the largest cluster. Right: Spanning 
probability. All systems were simulated for A = 1024 and £ =0.1. 

Strain rate dependence. While part of the y dependence 
was explored in Figure 2, we compare the transition curves for 
7= 10 3 and y = 10 6 in Figure S7. We observe a somewhat 
narrower transition for the lower strain rate, consistent with 
a sharp transition as y —> 0. (Not shown: Measures of the 
nonaffine parameter are consistent with this slight broadening 
at 7 = 10 5 ). Unfortunately, reliable N = 1024, y = 10 7 
results are beyond our current computational capacities, but 
results for smaller N point to a further narrowing. 

Other measures of rigidity. In rigidity percolation theory, an 
important measure is the spanning probability of the system, 
i.e. the probability of the existence of a system-size cluster. 
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FIG. S8. Probability of a spanning cluster for frictional (left) and 
frictionless (right) systems, same data set as Figure 4. 


This allows us to then pinpoint the unique transition point in 
the limit N —>• In Figure S8, we show the spanning proba- 
blities for three system sizes for both frictional (left) and fric¬ 
tionless (right) systems. For both, we observe a dramatic nar¬ 


rowing of the transition range with N. For both, we also con¬ 
firm that the transition point is in fact below the isostaticity 
point, at about dz = —0.15 for frictional systems and z = 3.9 
for frictionless ones. In Figure S7 (right), we also show the 
dependence of the spanning probability on y, again showing 
that higher strain rates slightly smooth out the transition. 
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